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Quantum algorithms for searching one or more marked items on a d-dimensional lattice provide an 
extension of Grover's search algorithm including a spatial component. We demonstrate that these 
lattice search algorithms can be viewed in terms of the level dynamics near an avoided crossing 
of a one-parameter family of quantum random walks. We give approximations for both the level- 
splitting at the avoided crossing and the effectively two-dimensional subspace of the full Hilbert 
space spanning the level crossing. This makes it possible to give the leading order behaviour for the 
search time and the localisation probability in the limit of large lattice size including the leading 

^^ order coefficients. For d = 2 and d = 3, these coefficients are calculated explicitly. Closed form 

^~H expressions are given for higher dimensions. 

O 

CN PACS numbers: 03.67.Ac, 03.65.Sq, 03.65.Aa 

'^ I. INTRODUCTION 

>--^ Quantum random walks as introduced by Aharonov, Davidovich and Zagury [1] in 1993 have gained considerable 

^^ attention over the last decade or so. It could be demonstrated that a quantum version of a classical random walk 

,—1 has transport properties which exhibit polynomial or even exponential speed-up compared to a classical random 

r^ walk, see [2-4] for overviews. This has recently lead to increased efforts for implementing quantum random walks 

C^ experimentally. In particular, one dimensional quantum walks have been realised using neutral atoms in optical 

-*— » lattices [5], with trapped ions [6-8], using coupled optical wave guides [9] and single photons [10]. 

C One of the most fascinating applications based on quantum random walk concepts are spatial quantum search 

p_^ algorithms. Like Grover's search algorithm [11, 12] for searching an unstructured data base, quantum walk search 

^^ algorithms can (usually) achieve a quadratic speed up compared to the corresponding classical search. The most 

1—^ prominent algorithms are the search on a hypercubc introduced by Shenvi, Kcmpe and Whaley [13] and the search 

on a d-dimcnsional lattice presented by Childs and Goldstone [14] in a continuous time version and by Ambainis, 

K^ Kempe and Rivosh (AKR) [15] as a discrete time algorithm. In particular, it has been pointed out that the search 

\^f-) in a two-dimensional lattice is critical with the number of search steps scaling like 0{\/N log N) whereas quantum 

f — ■ search algorithms on the hypercube as well as for lattices of dimension d> 3 scale like 0{-/N); here, N is the number 

^^ of vertices. Expressions for the leading order coefficients for the search time in the hypercube have been given in [16], 

^^ proposals for improving the efficiency of the hypercube search can be found in [17]. In [18], it has been shown that 

l/^ the number of steps for solving the search problem on a two-dimensional lattice can be decreased to 0{-\/W\ogN) by 

^^ modifying the quantum walk search algorithm, thus coming closer to the theoretical lower bound Vt{-\/N) [25] [19]. 

By changing the starting state, it has been demonstrated in [20] that quantum search algorithms on lattices can also 

be used in a sender-receiver configuration. Remarkably, this makes it possible to communicate across the lattice by 

sending information exclusively between two (or more) marked vertices where neither the sender nor the receivers 

need to know each others position. 

We will in the following focus on the AKR search algorithm on d-dimcnsional lattices. Extending the ideas and 
d techniques from [16], we will give improved estimates for the approximate eigenstate of the walk localised on the 
target vertex. This makes it possible to give closed form expressions for the leading order coefficients both for the 
search time and the search efficiency, that is, the overlap of the localised state with the target vertex. The paper 
is structured as follows: the search algorithm on the lattice and the basic vectors spanning the search space in the 
full Hilbert space are introduced in section II; this is followed by a detailed calculation of the normalisation constant 
for the approximate eigenvector localised on the target vertex in section III. In section IV, we analyse the spectral 
gap at the avoided crossing in an (approximately) invariant two-dimensional subspace. The results for leading order 
contributions to the localisation time and the amplitude at the target vertex are given in section V. 
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II. A QUANTUM SEARCH ON A d-DIMENSIONAL SQUARE LATTICE - INTRODUCING THE 

ALGORITHM 

Quantum search algorithms - Hke Grover's search - are usuaUy described in terms of a two unitary matrices: a 
quantum or wave propagator acting uniformly on the search space and a "marker" (or oracle) matrix which deviates 
from the identity only locally near the marked item. We will discuss here a lattice search algorithm of this type 
which was first introduce in [15]. We will limit the discussion to the search of one target vertex only. A generalisation 
to more than one marked item is straightforward and is discussed in [15, 20]. Wc will first introduce the propagator 
without a marked vertex leading to a quantum random walk on the lattice. The full quantum search algorithm and 
an analysis of the search mechanism will be given in Sec. II B. 



A. A quantum random walk on a d-dimensional lattice 

We consider a d-dimensional regular lattice with n vertices along each dimension, that is, N ^ nf^ vertices overall. 
The positions of the vertices in the lattice are defined by the set of vectors \x) = x = (a:;i,X2, . . . ,Xd) with integer 
coordinates Xi G {0, 1, . . . , n — 1}. Periodic boundary conditions will be assumed throughout. 

Following the spirit of quantum graph theory, we consider quantum waves propagating freely along 1-dimcnsional 
bonds connecting adjacent vertices [21, 22]. Waves undergo scattering at the vertices described by a 2fi x 2d scattering 
matrix a, where 2d is the number of bonds connected at each vertex. The scattering mechanism is the same at all 
vertices and the quantum dynamics on the lattice is in fact equivalent to a quantum random walk, see [23]. (The 
unitary scattering matrix a can then also be interpreted as a 'coin flip' matrix acting on an internal (spin-) degree 
of freedom of a quantum walker) . The total Hilbert space % of the problem is described as the tensor product of a 
vertex or position space \x) and a direction or coin space \v^) representing waves travelling in the negative/positive 
direction along the ith axis, i = 1, . . . , d. The dimension of T-L is thus 2dn'^ = 2dN . 

We chose a scattering matrix a which is unbiased with respect to the outgoing directions (except possibly for 
back-scattering). A natural choice is Grover's matrix [12] (which incidentally also features prominently in a quantum 
graph approach using Kirchhoff boundary conditions on the vertices, see [21, 22]). That is, we define [15] 

(T = 2|s)(s|-l2d, (1) 

where \s) is the uniform distribution in coin space, that is, 

^ i—l 

We have identical scattering processes at all vertices and we thus define a global coin flip matrix 



C = a®lN. (2) 

The wave propagation on the lattice is now given by a "shift" matrix S. Waves emanating from vertex \x) in the 
i th direction will reach vertex |x± e^) with |ei), the unit vector in the «^th direction. Writing \v^x) = \i ) (g) |a;), 
we find for the shift matrix S 

d 

S = y. /, {\i^x — Ci) {'i^x\ + \i^x + Ci) {i^x\) . (3) 

X i—l 

This is a natural choice for a moving shift since a walker leaving vertex x in direction i^ , enters vertex a; + e^ from 
direction i~ . 

The quantum random walk Uq is defined by first applying the global coin flip and then the moving shift 

Uo - SC. (4) 

In what follows, the eigenvectors and eigenvalues of Uq will be of importance which have been discussed in some 
detail in [15]. Using the tensor product, we write each eigenvector of Uq as a vector in coin space 
space vector IAj:) according to 



m| ) and a position 



^l) « \^k) 



where fc is a d-dimensional vector with components ki G {0, 1, . . . , n — 1} and c ~ 1, . . . 2d. The vector in position 

space can be factorised in the form IX^"^ = (2)^=1 IXfci)? where the \xki) sue obtained from the canonical basis vectors 
of position space using a Fourier transformation, that is, 

_. n— 1 

|;^fc.)= y^a'^'j'lj) where a = e2'^'/". (5) 

The thus obtained basis provides a convenient way to denote the eigenvectors in position space. 

For the quantum search, only IN — 1 of the 2dN eigenvectors of C/q are important, namely those having a coin 
space component not orthogonal to \s). These vectors are the 1-eigenvector |0o) = |s) ® P^g)' which is the uniform 
distribution, and two eigenvectors [15] 



«± 



■^.lA^i) (6) 

for each fc 7^ with complex conjugated eigenvalues e^'^^ and 

d 



2TTki 

cos tir = —, y cos 

d "^-^ n 



U-. ) in coin space become more and more cumbersome with increasing d and we 



The expressions for the vector 
will not attempt to give closed form expressions here. However, it has been shown in [15] that 

I "D - 75- <'' 

which is all we need in what follows. In addition, one has (d — 1) iV + 1 eigenvectors with eigenvalue 1 and {d— 1) N 
eigenvectors with eigenvalue —1 all with coin components perpendicular to \s). Note that the high degeneracies in 
the ±1 cigenspaces are due to the special choice of the coin flip matrices a and thus C. 

B. The quantum search algorithm 

1. The quantum search matrix 

We now assume that the lattice contains one marked vertex at (a yet unknown) position v; here, the marking is 
done by applying a different coin flip or scattering matrix a' at this target vertex |t;). The search algorithm is then 
defined as 

Ui ^ SC (8) 

C" = C - (cr - cr') ® \v) {v\ . (9) 

Note that, since \v) {v\ is a projection on the target vertex, the additional term changes at most Ad"^ matrix elements in 
C/q- As C/q is a 2dN x 2(iAf- matrix, Ui is identical to Uq up to a local perturbation and the quantum search algorithm 
can be regarded as a locally perturbed quantum random walk. 

Following AKR we define a' = — l2d which leads after some algebra to 

Ui^Uo{l-2\sv){sv\). (10) 

The perturbed quantum walk is shown in Fig. 1; starting from the uniformly distributed state |c^o) and applying 
Ui for t time steps with t = 0, 19, 38 and 57, one observes a localisation on the target vertex. Plotted here is the 
probability (wave amplitude square) on a logarithmic scale. 

Following the analysis developed for the hypercube in [16], we now define a one parameter family of unitary matrices 
Ux for the d-dimensional lattice: 

Ux^Uo + (c'^^ - 1) t/o \sv) {sv\ . (11) 
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Figure 1: Probability distribution of the quantum walk on a 31 x 31-lattice for a number of time steps t (on a logarithmic 
scale) . 

This makes it possible to extrapolate smoothly from the unperturbed walk (A = 0) to the AKR search algorithm 
(A = 1) and will be helpful for analysing the eigenvalues and eigenvectors of U\ near A = 1 in terms of the eigenbasis 
of the unperturbed quantum walk. 

Before doing so, we simplify the problem by reducing the size of the Hilbert space. We note that eigenvectors of 
Uq with eigenvalues ±1 (except the important vector \(j)o)) are orthogonal to \sv) and thus remain eigenvectors of U\ 
independent of the value of A; this part of H is thus irrelevant for the localisation effect. We define a reduced space H' 
by projecting out all ±l-eigenvectors of U except |0o). This subspace is a 2N — 1 dimensional Hilbert space spanned 
by |</.o)and|(/.i)for^^O. 

The scalar product ((/)- | sv) can be obtained using (7); thus an expansion of \sv) in terms of eigenvectors of Uq in 
the reduced space gives 

(12) 




Since \sv) is orthogonal to all eigenvectors of Uo not contained in H', this expansion also holds for the non-reduced 
space. 

2. Approximate eigenvectors of U\ 



In the following we will analyse the search algorithm in detail. Our approach is based on the work in [16] and makes 
it possible to go beyond the results found in [15] by giving explicit leading order coefficients for the search time. 

The overall strategy is best explained by considering briefly the spectrum of eigenphases {wj} of U\; in Fig. 2, 
the eigenphases are shown for an 11 x 11 lattice in the reduced space H'. The local perturbation leads to a set of 
avoided crossings in the spectrum, most notably at A = l,w = 0. This avoided crossing is formed by the state \(j)o) 
(which is a 1-eigenstate of U) and "another" state, which we will denote |z/a) in what follows. It will be argued that 
\vx) is localised in the vicinity of the target vertex and that the subspace {|<^o) i I^a)} is an approximate basis for the 
two-dimensional eigenspace spanned by the exact eigenstates at the crossing. The quantum search amounts then to 
a rotation in this two-dimensional subspace of "H' similar to the mechanism found in Grover's original search. It is 
important to note that unlike for Grover's algorithm, for the lattice search only approximate expressions for the state 
\v\=i) are known which will be given below; furthermore, the starting state |(/)o) has a small overlap with the rest of 
the spectrum at A = 1 leading to losses. 
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Figure 2: The eigenphases as functions of A for n = 11, d = 2. 



We start by considering the uniform distribution \(J)q) which is the remaining 1-eigenvector of Uq in "H'. One obtains 

1 



Ux |0o) = |0o> + (c'^" - 1) Uo \sv) 



N 



(13) 



that is, \(j>o) is an approximate eigenvector of U\ of order O (N^^'"^) . The remaining components, Uq \sv) , are locahsed 
on the 2d sites next to the target vertex. This can be seen from the definition of the shift operator S in Eq. (3) . 
To obtain Ivx), we start by expanding the vector in the basis of the unperturbed walk Uq, that is. 



a- 

k 



fe^O 

with a yet unknown set of coefficients {uq, a- }. Normafisation requires 



(14) 



|ao| 






We now make the ansatz 



Ux Wx) « e''?^^) \i.x) 



(15) 



(16) 



that is, we assume an approximate eigenvalue equation with a yet unknown eigenvalue e'^*-^'. This gives rise to a 
linear system of equations for the 2N — 1 unknown coefficients {ao,a-}. We are looking for a second approximate 
eigenvector which spans the space at the avoided crossing together with |(?!)o), that is we demand that \v\) and \(f>Q) 
are orthogonal which leads to oq = 0. 

A straight forward application of Ux on |i^a) results in 



UxWx)^J2{^P''' 

k^O 



'k,+-k^''' 



+ (e 



1 




sv I vx) 



with V, the position of the target vertex. In what follows, we set 

b := {sv I vx) 



(17) 



(18) 



and define the overall phase factor of \v\) such that h is real with < 6 < 1. Note that h determines the overlap of 
the approximate eigenvector with the target vertex. We will show in the next section that b is of order 0(l/\/\nN) 
for d = 2 and 0(1) for d > 3 at A = 1; thus \vi) is indeed localised at the target vertex. 
We now add and subtract the right hand side of Eq. (16) in (17) and obtain 






b{c' 



ittA 



N 



.igW 



a.\e 
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.i9(A) 



6(c 



l) a-'^Vs 



/2N 



b (c''^^ - l) a-^^e-'''R 



'2N 



(19) 



An approximate solution is obtained if each of the coefficients in front of the components 



6i ) is zero, that is, 



4 



b{e 



l) a-'=^e=t'*E 



/2N{&9W -e±'**E) 



(20) 



The resulting vector U\ \iy\) has a component in the \(J)q) direction of order 0{l/vN); thus, similar to Eq. (13), \iy\) 
remains in a subspace spanned by itself and \4>o). 

The coefficients ai enter Eq. (20) also through 6, that is, (20) represents a set of linear, coupled equations. Assuming 
6 7^ 0, equation (18) can be divided by b and the system of equations has a solution if and only if 



2iV 



E 

fe#0 



c'^aW 



cis'(^) 



-iOr 



(21) 



The avoided crossing occurs at A = 1 and c'^'^^ = 1, for which the right hand side of equation (21) can be obtained 
directly. 



-2 

2iV 



E(-l)-^(^-l)-l-^- 



(22) 



fe#0 



That is, no solution exists which maps a vector |j/i) exactly onto itself and \(f>o) ■ The error term 1/A^ is small compared 



to typical coupling terms between basis states |0*) (which are of the order l/VN) and the set of coefficients from 
Eq. (20) thus define an approximate solution of the eigenvalue equation (16). That is, the vector 

1) 



Wx) = 



b (e''^-^ 






27V 

— kv 



c^aW 



^k 



e 



iSr 



-■'Or 



(23) 



[26] fulfils to leading order the equation 



Ux Wx) - c's(^) I^^a) + 



5(< 



.ittA 



N 



(24) 



where 5(A) is defined implicitly by minimising the expression (21). An expansion of g(A) around A = 1 can be obtained 
successively, see [16]. We will restrict our attention here to the case A = 1 and thus |i^a=i) for which g(l) ~ 0. 

The only unknown quantity in determining |j^i) is the constant b which has dropped out of the equations. Since 
the vector |z^i) needs to be normalised, b turns out to be a normalisation constant. 



III. NORMALISATION OF li^i) 

Evaluating the normalisation constant b turns out to be the most laborious and technical part of the derivation. 
Readers more interested in the final results may want to proceed directly the result in equation (55). 

We will in the following restrict ourselves to the case A = 1 and thus e'^^^^ = 1 only. Demanding the vector in Eq. 
(23) to be normalised, one obtains 

"^ (25) 



\bf 2N^\i_^w, 



2 OAtA^.. ie,|2 



_2dy-^ 1 

" iVf-,d_Vf ,cos2^' ^ ' 

where the sum is taken over all vectors fc ^^ in a cf-dimensional cube with k,i G {0, 1, . . . , n — 1}. 

We start by rearranging the total sum into a summation over lower dimensional objects where only a limited 
number of entries of k are different from 0. That is, we consider the summation over all one dimensional edges, two 
dimensional faces, 3-dimensional cubes and so on, where the non-zero entries ki vary from 1 to n — 1. The edges 
are obtained by choosing rf — 1 entries of k equal 0, the faces have d — 2 entries of k equal and higher dimensions 
accordingly. This new arrangement results in 

±JAY '- (27) 

2rf . ^' 1 

a 



r^i 1 - eLi cos 



number of edges Jl 



contribution of one e( 

2d fd 



V 

,,,,=1 2 - E'=i cos : 



A^ \2j .^^^ '^- V^ ,f.,ns2ZL2i 

31,3 

number of faces 



contribution of one face 

+ . . . (28) 

^iV^(z). .^. . ,, _ V! , cos 2^ ^^^^ 



i=l ^"^ Jlj2,...Ji = l * Ei=lCOS „' 



i=l 



where li has been implicitly defined in (30) . Using the identity 1 — cos x ~ 2 sin | , we can simplify the sums and 
obtain 



(D" i: fe-'v - <-> 



Using Poisson's sunimation formula, one arrives at 



2 mi— — oo 

_i oo / i \~^ 

\x, J2 e^-™'^' E^i^'^ (32) 

mi = — DO \/=l / 

1 f^^"^ 



2 , dyi... I dyi 



°° °° g2iTi(mij/iH hniiyi) 



E ••• E %.. .. ' (33) 

E/=iSm yi 



?7ii — — OO mi— — OO 



where the last equahty is due to a reordering of terms and substituting yj — -Xj in ah integrals. 

The dominant contributions come from the terms tuj — for all j = 1, . . .i, that is, we write to leading order in 
the large n limit, 



1 



dyi 



dyi 



1 



'^J^ J^ Ei^iSin yi 

Using the symmetry of the sine squared function, the integrals (34) can be simplified to 

1 



/. = 2* 



dyi 



dyi 






(34) 



(35) 



We are only interested in leading order contributions and will thus use the integral (35) as an approximation for the 
full set of sums in Eq. (33). The expression (35) will form the starting point of our calculation of the leading order 
behaviour of \h\. 



A. Integration of I\ 



The first integration can be done explicitly and one obtains 

h = / dy—^^ = cot — = — + O - 
J JL. sm y 2n tt \n 



(36) 



B. Integration of lo 



To obtain a leading order estimate for 

/•f /-f 1 

I2 = 2 dx dy — 2 2~ 

J JL. J JL. sm a; + sm y 

we first note that the integrand is symmetric with respect to exchanging x and y and therefore 

1 



/2 = 4 / dx I dy^- 2 . 2 ■ 

J^ Jx sm x + sm y 

Now, the y integration is solved by observing that 

darctan(|^V2tan2a: + l) ^ 



Thus 



d^ sin .TV sin a; + 1 

arctan f^V2tan2a; 



• 2 ■ 2 ' 

sm X + sm y 



h = A dx 



2 2 
4 / dx- 



sinxv sin x + 1 
— arctan ( \/2 tan^ x + 1 



y=x 



sinxv sin^ x + 1 



Integrating by parts results in 



/2 =4 In ( J 2 + cot^ ^^ + cot ^^ I ( - - arctan J 2 tan^ ^^ + 1 



dx 



In ( \/2 + cot^ X + cot X j 



V2 + cot^ '. 



(37) 



(38) 



(39) 



(40) 



(41) 



(42) 



Using the substitution z ~ tana;, the resulting integration can finally be written in the forni[27] 

4 TT / 1 \ 

/2=7rlnn + 7rln--2i^--ln2 + C'^ , (43) 

TT 2 \n^ J 

where K « 0.916 is Catalan's constant. 

Using Eq. (30) , the result for the d — 2 dimensional lattice is 

1 4 v^ /2\ fn\i , , 

= -\nN+%{2-K)+'^\n%^o(^\. (45) 



C. Integration of I^ 

For the third integral /s, we will only evaluate the asymptotic limit for n — > oo. Starting with 

/•f /-f 1 

dx I dy dz^-^ ^^ „ , (46) 

J jL. J JL. sin a; + sin y + sin z 

we evidently have /3 > independent of n. An upper bound for I^ can be obtained using spherical coordinates, that 
is 

/3 < 4 / dr d^ de^-, —2 r^- (47) 

J jL^ Jo J^ sin X + sin y + sin z 

The integrand is bounded in the whole region of integration and in particular also in the limit r — > 0; we thus have 
h = Oil). 

Using numerical integration methods, we obtained 



lim I3 = 15.672 .... (48) 

n— ^00 

Substituting again into Eq. (30), we obtain for the normalisation of the d ^ 3 dimensional lattice: 



w^^iOGT'^ 



|2 

6 



- 3/3 = 6(1). (50) 

D. Integrations for d > 3 

We show in this section that 1/b^ — 0(1) for all d > 3; actual numerical values can be obtained by integrating 
expressions of the form (35) using numerical methods. From this, we can deduce that the search time T indeed scales 
like T ^ \/]V for d > 3 as will be shown in Sec. IV. 

We will proceed by induction: from (50) it is evident that 1/|6| = O (1) holds for d = 3. Thus, assuming that 
1/|6|2 = O (1) holds for some d > 3, it follows from Eq. (30) and li > 0, that Id = O (1). 

It remains to be shown that l/\b\ — O (1) also holds for d + 1. The starting point is again provided by Eq. (30), 
but this time for d + 1, that is, 

i,.2,.+i,(|:(^+')^,.-i,.,„). ,51) 

The sum of the first d terms adds to a leading order of 1/n, since li = O (1) for 3 < ^ < d and the lower order terms 
I2 = (In A^) and Ii = Q (n) have prefactors n^^'^ and n^'^, respectively. Therefore these contributions vanish for 



10 

It remains to be shown that I^+i = O (1) and this is done using Eq. (35), that is, 

Id+1^2^ f' dy,... f^ya+i—j r^ — . (52) 

As the squared sines are greater than zero, an upper bound is obtained by dropping the last sine term. Now, the j/^+i 
integration is performed, that is, 

= (' - s) '^ f^"' 

and this resuUs in J^+i = O (1). Thus, the normahsation constant b is of the order O (1). 
The overaU resuh for the normahsation constant is to leading order 



1 



\nN+ ^{2- K) + ^\n:^ for d = 2 
O (1) for d > 4 



^ I3 for d = 3 , (55) 



where K « 0.916 is Catalan's constant and we found I^ ~ 15.672 asymptotically, see (48). In particular, it has been 
shown that for i > 3, lim„_>.oo h converges. Thus, the leading order contribution to b can be calculated from equation 
(35) by replacing ^ with and using numerical integration methods. 

E. Localisation 

The normalisation parameter b is at the same time a measure for the localisation of |j^i) onto the state \sv) and 
thus on the target vertex v, see Eq. (18). Hence, the probability that the search will be localised at the target vertex 
\v) is proportional to 

|(z.l|s^;)|'-62. (56) 

(In fact, 6^ gives only a lower bound for the localisation probability as \s) fixes the coin space state in which the 
localised state is measured.). 

For estimating the overlap of the two-dimensional cigenspacc at the avoided crossing with the basis vector pair 
\(j)a), \i^i), we can resort to the results, Eqs. (13), (24). In particular, we found in (13), that the uniform state {(po) is 
mapped onto itself and a component in the direction Uo\sv). One finds by straightforward calculation 

(-1 I Uo I -) -—E { ,--J_,., + e-i.(i)-e-'^0 (5^) 

fe#0 

The overlap of |z^i) with Uq \sv), that is, the localisation on the nearest neighbours is thus of the same order as on 
the target state \sv) itself. To obtain the last equality e'^*^^^ = 1 has been used. 

We have shown that the vector |j/i) constructed in the preceding section is localised on the target vertex. We have 
furthermore shown that the vector space spanned by the orthonormal pair of approximate eigenvectors {\4>0i Wi)} has 
an order 0(1) or 0{l\/\nN) overlap with itself under the unitary map Ui. 

IV. QUANTUM SEARCH AT THE AVOIDED CROSSING 

The quantum search algorithm corresponds to a rotation in a lower dimensional sub-space spanned by the eigen- 
vectors at the avoided crossing at A = l,u; = 0. The search is initialised in the state \(j)o) which has an C(l) overlap 
with the exact eigenstates at the crossing, see (13). We have constructed a second approximate eigenvector at the 
crossing, namely \j^\=i) in (23), again with an 0(1) overlap with the eigenspace at the avoided crossing. We can thus 
restrict the analysis to studying the avoided crossing in the subspace of H spanned by the basis vectors {\4>o) , l^^i)}; 
this subspace is approximately invariant under the quantum walk Ui, as shown in Sec. II B, see also [16, 20]. 
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A. The size of the gap at the avoided crossing 

In a first step, we project out the two-dimensional submatrix related to the (approximately) invariant subspace, 
namely (C/^^^) = e^'^, where H is a, Hermitian 2x2 matrix. The entries of H can be determined by calculating the 
matrix elements of (t/^^^) explicitly. Using Eqs. (12), (13) and (24), one obtains to leading order 

(00 I C/i I 0o) = 1 - I (59) 

Wi \Ui\i^i) = l (60) 



for the diagonal elements and 



6o I C/i I ^i) = -^ (61) 



(t^i \Ui\4>o)^-^ {i^i I Uo I sv) (62) 



^(-^1 



for the off-diagonal entries, where Eq. (58) has been used for the last equality. Identifying the basis vectors {{ipa) , Wi)} 
with {|1) , 1 2)}, we obtain to leading order 



r2x2 /'l -eA •-- 26 



^1=1. 1 ^i*^ '=^ (64) 



^ - [t 7) • (65) 



which results in a Hamiltonian 



We note that the coupling term e scales like 1/\N and depends on the normalisation constant 6. 
Eigenvectors and eigenvalues of H are easily calculated, that is, 

eigenvector eigenvalue 

1 A' 



1"^^ = 73 \^) -' (66) 



This means in particular that the gap at the avoided crossing behaves like 



The time of search is directly related to the gap A; the search starts in the vector |(/)o) = l/\/2(|ui) + |m2)) which 
after successive iterations becomes 

(C/f 2)* |</-o) = c-^* |</-o) = c-^*-^ {\u{) + \u^)) (67) 

= -^(e-*|«i)+e--*|«2)). (68) 
Thus, for a time ?" •= ^ — s, one has e*"^* — ±i and 

(L/r')*|0o) = ^(K)-K)) = -|^i), (69) 

which is the localised state. Note, that the search time T is inversely proportional to the coupling parameter e in Eq. 
(65) and thus to the spectral gap at the avoided crossing. The quantum search algorithm succeeds after 0{VN/b) 
steps and thus faster as any classical search. 
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Figure 3: The localisation time T after which the walk is localised at the marked vertex for several A''. Numerical result 
compared to analytical for d = 2 and 3. 



RESULTS 



Putting everything together, we have now identified the search time in leading order as 

TT TT\/N 



T 



2e 



46 



(70) 



It depends on the normalisation constant 6 which has been obtained in (55). Overall, the number of time steps T 
can now be given as the integer closest to 



^^/fln7V+A(2-i^) + lln^ ford = 2 

for d = 3 
Q[^JN\ for d > 4. 



(71) 



Note, that the bound for d > 3 is tight, because we know from [11, 19], that the search can not be faster than 
VW. In Fig. 3, we compare the analytical results for the localisation time with numerical simulations. The theoretical 
results agree very well with the outcome of the numerical simulations. The general behaviour suggests that the walk 
for a fixed number of vertices is the faster the higher the dimension. 

The probability to find the search at the target vertex after T steps is equal to 6^ (18). From Eq. (55), we find 



ii 



lnW^ 



■In: 



■ln2+0(i) 



e 



/InN 



i-. 



.e(i) 



9(1) 



for d == 2 
for d = 3 
for d > 4 



(72) 



Note, that the upper bound for b in the case d > 4 is a consequence of (70) and the lower bound for the search time 
is provided by [19]. 

To leading order, the localisation amplitude at the target verte x is n ot A^-dependent for d > 3; the d = 2 case 
is exceptional, as the projection of |j^i) onto \v) decreases like l/VlniV. Note, however, that the state |j^i) is still 
localised on the marked vertex for d ~ 2 since the average amplitude of an eigenstate of the random walk Uq on an 
arbitrary vertex decreases like -4=. Amplitude amplification methods can be used to increase the success probability 

to a constant in N by repeating the search algorithm O ( y/hiN) times [24]. 

In Fig. 4, we compare the numerical results for the probability to find the walk at the target vertex with the analytical 
results for 6^. The numerical results are obtained from running the quantum search algorithm and determining the 
maximum of the probability at the target vertex. It is noted that for d = 3, the asymptotic value given in Eq. (72) is 
not reached before N « 10^, (see inset of Fig. 4). This corresponds to a 3d lattice of side- length n = 100 - far beyond 
what can be simulated on a computer. However, a direct evaluation of b^ in terms of the sums in Eq. (30) confirms 
the basic mechanism behind our approach. 
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Figure 4: Maximal probability at the target vertex as a function of A*': circles {d = 2) and squares (d — 3) are obtained directly 
from the quantum search; the solid and dashed lines represent Eq. (72). The dashed-dotted line gives the result for b^ obtained 
directly from Eq. (30) for d — 3. Inset: Direct evaluation of Eq. (30) compared with asymptotic result Eq. (72) for d — 3 over 
large N interval. 

VI. CONCLUSIONS 

We presented a detailed analysis of the lattice search and shovi? in particular that the search mechanism is based on 
a rotation from a uniformly extended state to a localised state coupled at an avoided crossing point in the spectrum 
of the operator U\ . The search time and localisation strength can be obtained from an effectively two-dimensional 
model similar to the analysis in [16]. We have thus independently verified the scaling behaviour of the search time and 
localisation probability as stated in [15] and can now also give explicit expressions for the leading order coefficients; 
the results have been verified by numerical simulations. 

A generalisation to a search with more than one marked item is straightforward and leads to reduced models as 
presented in [20]. We also give a detailed derivation of how to obtain the localised state \i^\) asymptotically. This 
is important for the wave communication protocols proposed in [20]; here, the search algorithm is used to search for 
more than one marked vertex, to send signals through the lattice from a sender to one or more receiver points and to 
construct new search algorithms that do not rely on knowing the number of target vertices. 
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